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SUMMARY 


Details  are  given  of  calculations  of  the  flow  with  detached  bow 
wave  past  a  doubly  symmetrical,  double-wedge  profile  at  zero  angle  of 
attack.  The  results  of  these  calculations  have  been  described  pre¬ 
viously  in  MCA  TN  2339-  The  mixed  flow  over  the  front  of  the  profile 
is  determined  by  the  solution  of  a  boundary- value  problem  for  the 
transonic  small -disturbance  equation  in  the  hodograph  plane.  The  solu¬ 
tion  is  obtained  by  the  use  of  finite-difference  equations  and  relax¬ 
ation  techniques.  The  methods  follow  established  lines  except  for  the 
somewhat  novel  treatment  of  the  boundary  conditions  along  the  shock 
polar  and  sonic  line.  The  purely  supersonic  flow  over  the  rear  of  the 
profile  is  determined  by  means  of  the  method  of  characteristics  special 
ized  to  the  transonic  small-disturbance  theory. 


INTRODUCTION 

The  present  report  is  the  second  of  two  NACA  publications  on  a 
theoretical  study  of  the  transonic  flow  past  a  wedge  profile  with 
detached  bow  wave.  The  specific  problem  under  consideration  is  that 
of  a  thin,  finite'  wedge  at  zero  angle  of  attack,  with  particular  appli¬ 
cation  to  the  case  of  the  doubly  symmetrical,  double-wedge  profile. 

The  first  report  (reference  l)  contains  a  nonmathematical  outline  of 
the  general  analytical  method  and  a  discussion  of  the  final  calculated 
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results.  The  present  report  supplements  the  previous  work  by  supplying 
the  details  of  the  analysis.1 

The  plan  of  the  report  is  briefly  as  follows:  In  the  first  section, 
the  basic  problem  of  the  finite  wedge  with  detached  bow  wave  is  stated 
as  a  boundary-value  problem  for  the  transonic  small -disturbance  equation 
in  the  hodograph  plane.  Except  for  the  introduction  of  a  boundary  con¬ 
dition  along  the  sonic  line  to  replace  conditions  previously  prescribed 
In  the  supersonic  portion  of  the  hodograph,  this  material  follows  th§ 
lines  established  by  Frankl  (reference  k)  and  Guderley  (reference  5). 

It  is  recounted  here  primarily  for  the  sake  of  completeness .  In  the 
second  section  of  the  report,  the  boundary-value  problem  in  the  hodo¬ 
graph  is  reduced  to  a  system  of  finite-difference  equations  the  solution 
of  which  is  then  obtained  by  relaxation  techniques.  This  portion  of  the 
work,  which  constitutes  the  main  contribution  of  the  present  report,  is 
discussed  in  some  detail,  since  it  is  anticipated  that  the  methods  and 
equations  which  are  presented  will  be  useful  in  the  solution  of  other 
problems  involving  detached  shock  waves.  The  third  section  of  the 
report  describes  the  transformation  of  the  hodograph  solution  for  the 
finite  wedge  back  into  the  physical  plane.  The  fourth  section  is  con¬ 
cerned  with  the  characteristics  construction  used  to  obtain  the  purely 
supersonic  flow  over  the  rear  of  the  double-wedge  profile,  and  the  final 
section  contains  a  few  remarks  on  the  accuracy  of  the  solution.  It 
will  be  presumed  that  the  reader  is  familiar  with  the  general  discussion 
of  the  investigation  given  in  reference  1  or  with  the  description  of  the 
finite-wedge  problem  available  in  reference  5>  6,  or  7*  2 


!As  explained  in  the  introduction  to  reference  1,  it  was  originally 
intended  that  a  third  report  be  included  under  the  present  general 
title.  This  report,  written  by  Arthur  E.  Bryson,  Jr.,  was  to  describe 
an  experimental  study  of  wedge  profiles  conducted  at  the  California 
Institute  of  Technology.  Because  of  a  broadening  of  the  experimental 
work:  to  include  biconvex  profiles  and  subsonic  as  well  as  supersonic 
flight  speeds,  it  now  seems  advisable  to  depart  from  the  original 
plan.  BrysonT  s  report  is  therefore  being  issued  by  the  NACA  under  a 
less  restrictive  title  (reference  2).  Certain  of  the  results  of  both 
the  theoretical  and  experimental  studies  have  been  given  in  prelimini¬ 
nary  form  by  Liepmann  and  Bryson  in  reference  3* 

2The  authors  are  indebted  to  William  A.  Mersman  of  the  Ames  Laboratory 
for  suggestions  leading  to  certain  of  the  mathematical  procedures  used 
in  the  analysis. 
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NOTATION 
Primary  Symbols 

a*  critical  speed  (i.e.,  speed  at  which  the  speed  of  flow  and 

the  speed  of  sound  are  equal) 

c  airfoil  chord 

e >f >h,k  length  of  irregular  lattice  intervals  relative  to  that  of 

the  basic  interval 

F(\)jG(A,)  functions  defined  along  sonic  line  in  hodograph 
(See  equation  (6).) 

^w  integral  defined  by  equation  (39) 

<Jj_  (i  =  ij2,  3)  component  integrals 

(See  equation  (33)  et  seq.) 

K,L  functions  of  e,  f,  and  r)0 

(See  equation  (27).) 

constants 

(See  equations  (7)  and  (15).) 

M  Mach  number 

S(tJ)  function  defined  by  equation  (24) 

t  airfoil  thickness 

V  local  speed  of  flow 

x,y  Cartesian  coordinates 

Y  ordinate  function 

(See  page  31.) 

3  absolute  value  of  “rj  at  left-hand  limit  of  lattice 

J  ratio  of  specific  heats  (1.4  for  air) 

A  basic  lattice  interval 

T)  ,  hodograph  variable  defined  by  equation  (l) 
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local  inclination  of  flow 
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■* 


0' 


®u>©d 


X 


Xf 


£ 


£0 


P 


* 

+E 


variable  of  integration 
(See  equation  (8).) 

semiangle  of  wedge 

ordinates  of  upgoing  and  downgoing  characteristics  at  *tj  =  0 
(See  page  37.) 

hodograph  variable 
(See  equation  (6).) 

variable  of  integration 
(See  equation  (6).) 

speed  parameter 

(See  equation  (49).) 

transonic  similarity  parameter 
(See  equation  (17).) 

fluid  density 

stream  function 

value  of  at  the  point  E 
(See  equation  (l6).) 


Subscripts 


o  conditions  in  free  stream 

o,i ,  2, etc.  value  at  a  prescribed  lattice  point 

*  conditions  at  critical  speed 


Superscript 


quantity  in  normalized  form 
(See  equation  (9)0 


A* 
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STATEMENT  OF  BOUNDARY -VALUE  PROBLEM  IN 
HODOGRAPH  PLANE 

Complete  Holograph 


A  qualitative  representation  of  the  flow  over  a  double-wedge 
profile  with  detached  bow  wave  is  given  in  the  first  of  the  accompanying 
sketches,  which  differs  only  in  detail  from  a  similar  sketch  presented 
in  reference  1.  A  description  of  this  flow  field  is  to  be  found  in  the 


Vo/a,  >  / 


-  Shock  waves 

- Sonic  line  (V/a,=l) 

— » -  Streamline 

- Expansion  \  Mach 

- —  Compression  /  lines 


earlier  paper.  The  second  sketch  shows  the  corresponding  picture  for 
the  flow  about  the  front  half  of  the  profile  in  the  holograph  plane  - 


that  is,  in  the  plane  V/a*,  0  where  V  is  the  local  speed  of  flow, 
a*  is  the  critical  speed,  and  6  is  the  local  inclination  of  flow 
relative  to  the  x  axis. 


The  picture  in  the  upper  half  of  the  holograph  plane  can  be 
described  briefly  as  follows:  The  part  of  the  shock  wave  which  borders 
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upon  the  subsonic  region  in  the  physical  plane  appears  in  the  hodograph 
as  the  subsonic  portion  AE  of  a  shock  polar.  The  shape  and  position 
of  the  shock  polar  are  determined  by  the  dimensionless  free-stream 
velocity  V0/a*  (or,  what  is  equivalent,  by  the  free-stream  Mach  number 
M0)  and  by  the  ratio  of  specific  heats  7. 3  The  portion  of  the  central 
streamline  from  the  normal  part  of  the  shock  wave  to  the  stagnation 
point  at  the  nose  of  the  wedge  maps  into  the  portion  AO  of  the  hori¬ 
zontal  axis  in  the  hodograph.  The  image  of  the  wedge  itself  is  given 
by  a  radial  line  inclined  at  the  wedge  angle  0W  and  extending  from 
the  origin  0  to  the  point  B  on  the  critical  circle  (v/a*  =  l) . 

The  shoulder  of  the  wedge,  which  produces  an  expansion  fan  of  a  locally 
Prandtl-Meyer  type  in  the  physical  plane,  appears  in  the  hodograph  as  a 
portion  of  the  downgoing  characteristic  (epicycloid)  starting  at  B. 

The  last  Mach  line  from  the  shoulder  to  the  sonic  line  (termed  the 
separating  Mach  line  in  reference  l)  appears  as  part  of  the  upgoing 
characteristic  which  begins  at  the  intersection  E  of  the  shock  polar 
and  critical  circle.  Point  G,  the  point  of  intersection  of  the  epi¬ 
cycloids  from  B  and  E,  fixes  the  extent  of  the  downgoing  character¬ 
istic  which  must  be  considered  in  determining  the  solution  in  the 
hodograph.  A  typical  streamline  in  the  hodograph  plane  is  shown  by  the 
line  FCD. 

To  obtain  a  solution  for  the  detached-wave  problem  in  the  hodograph 
a  boundary- value  problem  for  the  differential  equations  of  gas  dynamics 
must  be  solved  within  the  region  AOBGEA.  If  the  stream  function  + 
is  taken  as  the  unknown,  the  pertinent  boundary  conditions  are  as 
follows: 

1.  The  value  of  is  constant  along  the  basic  streamline  AOBG. 

2.  The  streamlines  (i.e.,  the  lines  of  constant  t)  leave  the 
shock  polar  with  a  direction  which  is  a  known  function  of  location  on 
the  polar. 

3*  The  increment  in  ^  over  the  portion  AE  of  the  shock  polar 
has  a  prescribed  value  different  from  zero. 

The  reason  for  the  first  condition  is  obvious.  The  second  condition  is 
a  consequence  of  the  requirement  that,  at  every  point  on  the  shock  polar 
the  direction  of  the  shock  wave  as  computed  from  the  solution  for  \|r 
must  be  compatible  with  the  direction  given  by  the  equations  for  an 
oblique  shock  wave  in  a  uniform  stream.  The  third  condition  prevents  ^ 
from  being  simply  a  constant  throughout  the  hodograph  and,  in  effect, 
fixes  the  scale  of  the  flow  field  in  the  physical  plane.  It  will  be 
noted  that  no  condition  is  prescribed  along  the  boundary  EG  in  the 
hodograph.  Frankl  has  proved  (reference  4)  that  the  solution  determined 
by  the  foregoing  boundary  conditions  is  unique. 


3 The  equations  which  are  pertinent  here  can  be  found  in  the  work  of 
Frankl  (reference  4). 
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Specialization  to  Small  Disturbances 


Original  boundary- value  problem,  -  As  has  been  shown  by  Guderley 
(reference  5)^  the  equations  of  the  boundary-value  problem  in  the  hodo¬ 
graph  are  considerably  simplified  when  restriction  is  made  to  the 
neighborhood  of  the  critical  speed*  To  this  end,  the  quantity  r\  is 
introduced  according  to  the  relation 

,  V-a* 

n  =  (7+i)1/3  —  (i) 

a*  f 

t 

and  T)  and  the  stream  angle  Q  are  assumed  sufficiently  small  that 
only  their  lowest  powers  need  be  retained  in  the  analysis.  This  means, 
in  effect,  that  the  right-hand  portion  of  the  previous  hodograph 
(including  the  shock  polar  itself)  is  made  to  shrink  down  to  the  vicin¬ 
ity  of  the  point  H,  which  defines  the  intersection  of  the  critical 
circle  and  the  horizontal  axis.  ' 

When  the  foregoing  procedure  is  carried  out  and  the  limiting 
process  is  counteracted  by  a  suitable  enlargement  of  scale,  the  situ¬ 
ation  in  the  small-disturbance  hodograph  (i.e.,  in  the  tj,0  plane) 
appears  as  in  the  adjacent  sketch.  Here  the  critical  speed  corresponds 
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to  the  vertical  axis  tj  =  0.  The  equation  for  the  upper  half  of  the 
shock  polar  in  the  simplified  hodograph  has  the  form 

»  •  <%-1>  (2) 
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•where  T|0  ■  is  the  value  of  tj  corresponding  to  the  free-stream 
velocity  VQ.4  By  virtue  of  the  limiting  process,  the  stagnation  point 
at  0  has  moved,  in  the  present  system  of  axes,  infinitely  far  to  the 
left.  As  a  result,  the  part  AO  of  the  horizontal  axis  (6  =  0)  extends 
now  from  f\  =  -T)0  to  T)  =  -<».  The  image  OB  of  the  wedge  is  similarly 
represented  hy  the  horizontal  line  9  =  0-y,  T)  =  s0.  The  characteristics, 
which  complete  the  boundaries  of  the  field,  have  the  simple  form 


9  =  const. 


±  -  n3/2 

3 


(3) 


On  the  basis  of  the  usual  assumptions  regarding  flow  near  the 
critical  speed,  the  differential  equation  for  reduces,  in  the  pres¬ 
ent  simplified  hodograph,  to  the  form 


+TJTI  -  ^ee  58  0 

This  is  the  mixed  elliptic-hyperbolic  equation  studied  by  Tricomi  in 
reference  8.  The  boundary  conditions  along  the  central  streamline 
require  that  is  constant'  -  say  0  -  on  AO  and  OBG  and  that 

i|r->0  as  T]->(-co)  for  0^0^  0W.  On  the  upper  half  of  the  shock  polar, 
the  boundary  conditions  require  that  the  lines  of  constant  must  have 
the  slope 

d0  =  _  t)q+7ti  Ao+Tj  (5a) 

dTj  3ti0+5t1  *  2 

On  a  line  of  constant  4”,  d@/dT|  can  be  replaced  by  /^0 ,  so 

the  foregoing  condition  can  also  be  written 


Min  =  0 

3ti0+5t)  v  2  9 


(5’d) 


The  final  boundary  condition  requires  that  must  have  some  given 
value  ihr(A>)  at  the  point  E.  Since  the  coordinates  of  the  flow  field 
uil i  ultimately  be  expressed  in  terms  of  a  characteristic  dimension  of 
the  wedge,  the  actual  value  assigned  to  is  purely  a  matter  of  con¬ 

venience.  As  before,  no  boundary  condition  is  specified  along  the 
characteristic  EG. 


4 The  derivation  of  this  and  the  other  equations  for  the  simplified 
hodograph  is  given  by  Guderley  in  reference  5* 
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Elimination  of  the  supersonic  region.-  The  foregoing  is  the 
boundary-value  problem  for  the  finite  wedge  as  formulated  by  Guderley. 

It  was  the  original  intention  in  the  present  work  to  obtain  a  numerical 
solution  of  this  problem  on  the  basis  of  the  boundaries  and  boundary 
conditions  which  have  been  described.  Efforts  in  this  direction  failed, 
however,  because  of  difficulties  in  obtaining  convergence  of  the  relax¬ 
ation  process  in  the  supersonic  portion  BGE  of  the  holograph.5 
Similar  difficulties  have  been  reported  in  references  9  and-  10  with 
regard  to  relaxation  calculations  of  the  transonic  flow  through  a 
converging-diverging  nozzle.  The  reasons  for  the  difficulty  in  the 
present  case  are  not  apparent.  Fundamental  .questions  would  appear  to 
be  involved  concerning  the  stability  and  convergence  of  the  finite- 
difference  scheme  for  the  Tricomi  equation  in  the  hyperbolic  domain.  A 
study  of  these  matters,  similar  perhaps  to  that  reported  for  the  wave 
equation  in  reference  11,  may  be  a  prerequisite  to  numerical  solutions 
of  mixed-flow  problems  in  the  general  case.  In  the  present  example, 
however,  the  difficulty  can  be  circumvented  by  modifying  the  boundary- 
value  problem  so  as  to  eliminate  the  supersonic  region  from  explicit 
consideration. 


The  elimination  of  the  supersonic  region  depends  on  a  formula 
given  by  Tricomi  (reference  8,  equation  (2.19))  which  relates  the 
behavior  of  \(r  on  the  vertical  axis  to  its  behavior  on  a  character¬ 
istic.  In  the  present  case,  in  which  ^  is  identically  zero  on  the 
characteristic  BG,  this  formula  reduces  to  an  integral  relation 
between  i|f  and  ^  at  points  on  the  sonic  line.  This  relation  has 
the  form 


F(\) 


J  (x-x')1/3 


where 


X  =  Q-y  -  Q 

f(x)  =  t(o,e) 

G(x)  =  ^(0,0) 

and  X1  is  a  variable  of 
integration.  (See  sketch. 


G 


5Several  procedures  were  tried  in  the  supersonic  region,  using  both  a 
square  lattice  and  a  lattice  following  the  characteristics.  All  were 
unsuccessful. 
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The  numerical  constant  kx  is  given  ty 


.  t!±  r»(l/3)  (7a) 

where  T(l/3)  is  the  gamma  function  of  the  argument  1/3.  For  the 
accuracy  required  in  the  later  numerical  work,  the  value  of  this  con¬ 
stant  may  be  taken  as 

k-L  =  1.013  (7b) 

Satisfying  equation  (6)  everywhere  on  the  sonic  line  from  B  to  E  is 
completely  equivalent  to  satisfying  the  condition  ♦  =  0  on  the  down¬ 
going  characteristic  from  B  to  G. 

For  the  present  application,  it  is  convenient  to  invert  equa¬ 
tion  (6),  which  can  be  done  by  means  of  Abel’s  formula.  This  gives 
(see  reference  12,  p.  229) 


0(X)  =  fX  F<X'>  dX' 

2itk1  ax  J  (X-X')2/3 

The  differentiation  indicated  on  the  right  is  readily  accomplished  by 
first  transforming  the  integral  to  one  with  fixed  limits  by  means  of 
the  substitution  Xf  =  t\.  The  result  is,  after  reverting  to  the  orig¬ 
inal  notation. 


G(X) 


VT  1  r  1 

2rtkx  X  L  3 


fx  F<v>  ax'  +  rx  bLMn.  ax' 
J  (X-X*)2/3  J  (X-X*)2/3 


Transforming  the  first  integral  through  integration  by  parts  and  noting 
that  F(o)  =  0,  one  obtains  finally 


G(X) 


71  fx  _£xM_  ax. 

2rtka  J0  (X-X«)2/3 
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This  can  be  written  in  the  rj,0  notation  as  * 


^(0,8)  + 


*6(0,9') 

(e»-0)2/3 


dS» 


=  0 


(8) 


where  G'  =  G^-X'  denotes  the  variable  of  integration.  As  with 
equation  (6),  satisfaction  of  equation  (8)  everywhere  on  the  sonic 
line  between  B  and  E  insures  that  ^  is  zero  everywhere  on  the 
characteristic  from  B  to  G.  By  regarding  equation  (8)  as  a 
boundary  condition  along  BE,  the  region  of  solution  of  the  partial 
differential  equation  (4)  can  be  confined  to  the  purely  subsonic 
portion  of  the  holograph  (t)^0).  Relaxation  methods  can  be  used  to 
solve  the  resulting  elliptic  problem  without  essential  difficulty. 

Equations  in  normalized  form.-  To  carry  out  the  numerical 
calculations,  it  is  convenient  to  normalize  the  equations  of  the 
boundary- value  problem  by  means  of  the  transformation 


»o3/2 


(9) 


This  is  equivalent  to  introducing  the  rules  for  transonic  similarity. 
(See,  for  example,  references  5,  13#  and  lit  - )  The  particular  form 
of  transformation  chosen  here  has  the  advantage  for  the  present  work 
of  providing  a  unique  shock  polar  with  conveniently  located  hori¬ 
zontal  and  vertical  intercepts. 

With  the  foregoing  substitution,  the  differential  equation  (b) 
takes  on  the  following  form  in  the  %(?  plane: 


r\  T) 


=  0 


(10) 
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On  the  sonic  line  BE: 


uo,s)  +  k2  re 

n  J*.  (e'-e)2/3 


iSe<e. 


where 


k2  = 


3 1/2 
22/3itk1 


4.  At  the  point  E: 


=  i|rE  (arbitrary  constant  ^0)  for  1^0,  $  =  1  (16) 

It  is  apparent  from  the  preceding  equations  that  a  solution  of  the 
problem  will  depend  on  only  the  single  ^parameter  which  defines  the 

position  of  the  upper  boundary  in  the  tJ,0  plane.  This  parameter  is 
directly  related  to  the  more  familiar  transonic  similarity  parameter 

6  « _ 2^ _  (17) 

0  [(7+l)(t/c)]2/3 

which  was  used  as  a  basis  for  the  presentation  of  the  results  in 
reference  1.  (in  this  equation,  t/c  is  the  thickness  ratio  of  the 
double-wedge  profile.  It  is  equal,  to  the  present  order  of  accuracy, 
to  the  wedge  angle  0W.)  The  relation  between  |0  and  0W,  easily 
derived  from  equations  (l)  and  (9)  of  the  present  paper  and  equa¬ 
tion  (2l)  of  reference  1,  is  given  by 


,  =  2V3 
0  V2/3 


SOLUTION  OF  BOUNDARY-VALUE  PROBLEM 
IN  HODOGRAPH  PLANE 


The  solution  of  the  boundary-value  problem  in  the  T},0  plane  is 
obtained  'in  two  steps,  according  to  established  procedures  for  the 
numerical  treatment  of  partial  differential  equations.  (For  intro¬ 
ductory  articles,  see  references  15,  16,  and  17 .  For  an  extended 
discussion,  see  reference  18.)  In  the  first  step,  the  domain  Under 
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consideration  is  covered  by  a  square  lattice,  and  a  finite-difference 
approximation  to  the  differential  equation  or  boundary  condition  is 
written  for  each  lattice  point.  The  boundary- value  problem  for  the  par¬ 
tial  differential  equation  is  thus  reduced  to  a  problem  in  solving  a 
large  number  of  simultaneous  algebraic  equations.  Solution  of  the 
latter  problem  by  relaxation  methods  is  the  second  step. 


Reduction  to  Finite-Difference  Equations 


The  arrangement  of  a  typical  finite-difference  lattice  in  the  ^,0 
plane  is  shown  in  figure  1.  The  basic  lattice  interval,  which  is  the 
same  in  both  directions,  is  denoted  by  A.  Adjacent  to  the  shock  polar, 
the  interval  is  adjusted  so  that  the  terminal  lattice  points  lie  on  the 
polar  itself.  For  purposes  of  formulating  the  finite-difference  equa¬ 
tions,  the  lattice  points  are  conveniently  grouped  into  five  categories 
as  follows  (typical  points  in  each  category  are  indicated  in  the  figure): 

a.  Regular  internal  points 

b .  Points  far  to  the  left 

c.  Points  adjacent  to  the  shock  polar 

d.  Points  on  the  shock  polar 

e.  Points  on  the  sonic  line 

The  form  of  the  finite-difference  equation  pertinent  to  each 
category  will  be  developed  in  the  subsequent  paragraphs.  The  methods 
employed  are  standard,  except  for  the  somewhat  novel  treatment  of  the 
boundary  conditions  along  the  shock  polar  and  sonic  line. 

Regular  internal  points.-  The  category  of  regular  internal  points 
comprises  all  points  interior  to  the  boundaries  but  not  immediately 
adjacent  to  the  shock  polar.  The  situation  in  the  vicinity  of  such  a 

point  is  as  shown  in  the  sketch.  The  differ¬ 
ence  equation  which  applies  here  is  obtained 
by  suitable  approximation  to  the  differential 
equation  (10). 

If  it  is  assumed  that  the  unknown 

/v 

Z  function  i|r  =  i|f(r),0)  may  be  expanded  locally 
in  the  form  of  a  Taylor’s  series,  the  values 
of  at  2  and  4  may  be  written  (see  sketch) 
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^2 

=  to  +  At- 

+  — 

+  o(a4) 

n  0 

21 

TIT)  0  31 

T)T)T)  0 

*4 

=  -  At 

0  * 

t  I 

+  o(a4) 

T)  O 

21 

nn  p  31 

mjp 

Addition  of  these  equations  and  solution  for  \jr^i  gives 

Vl|o 

♦—I  =  \  (^4-2^0+^)  +  °(a2)  (19a) 

in|o  a2 

which  is  a  well-known  difference  egression  for  the  second  derivative. 
The  corresponding  derivative  in  the  vertical  direction  is  similarly- 
represented  by 

tv~  =  h  (VZV*i)  +  °(a2)  (!9b) 

00o  A 

Substituting  these  expressions  into  the  differential  equation  (lO)  and 
neglecting  the  terms  of  0(A2)  then  gives  for  the  finite-difference 
equation  at  a  regular  internal  point 

t-a  +  +4  “  ^»0(1lri+1lrs)  "  2(1-2?) oH0  =  0  (20) 

where  Tj0  denotes  the  horizontal  coordinate  of  the  point  in  question. 
The  difference  equation  for  internal  points  is  thus  the  same  for  points 
on  a  given  column  but  differs  from  one  column  to  the  next.  For  a  point 
adjacent  to  the  upper  boundary  OB  —  as,  for  example,  the  point  a1  in 
figure  1  —  the  value  of  ^  must  be  set  equal  to  zero  in  accord  with 
the  boundary  condition  (12)  .  Similar  considerations  hold  for  points 
adjacent  to  the  lower  boundary  OA. 

Equation  (20)  represents  the  simplest  possible  finite -difference 
approximation  to  the  differential  equation  (lO)_.  As  is  apparent  from 
the  derivation,  the  error  involved  is  of  0(A2) .  Consideration  has 
been  given  to  improving  the  approximation  by  including  additional  lat¬ 
tice  points  in  the  finite-difference  equation  or  by  incorporating 
higher-order  difference  corrections  in  the  later  relaxation  work  (cf., 
references  19,  20,  and  21) .  Because  of  the  complicated  nature  of  the 
boundary  conditions  along  the  shock  polar  and  sonic  line,  however, 
consistent  application  of  these  procedures  did  not  appear  feasible. 

The  requisite  accuracy  in  the  present  work  has  therefore  been  achieved 
by  suitable  decrease  in  the  mesh  interval  A  in  those  regions  in  which 
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the  function  ^  varies  most  rapidly.  This  procedure  has  the  secondary 
advantage  of  providing  closely  spaced  values  of  the  derivatives  which 
are  required  for  the  later  transformation  to  the  physical  plane. 

Points  far  to  the  left.-  In  order  to  carry  through  the  numerical 
analysis,  it  is  necessary  that  the  finite-difference  lattice  be  termi¬ 
nated  at  some  distance  to  the  left  in  the  hodograph.  This  can  be  done 
with  the  aid  of  an  asymptotic  solution  valid  for  large  negative  values 
Of  Tf. 


By  separation  of  variables,  it  can  be  shown  that  the  general 
solution  of  the  differential  equation  (lO)  in  the  region  -1, 

0  =  $w>  subject  to  the  boundary  conditions  (ll),  (12),  and  (13)>  is 


00  _  /v 

4?  (-2T))3/2" 

3@w 

n=i 

where  Kx /3  is  the  modified  Bessel  function  of  the  second  kind  of  order 
1/3  (notation  of  reference  22)  and  An  is  an  appropriate  constant.  At 
sufficiently  large  negative  values  of  rf  the  first  term  of  the  series 
will  predominate,  and  the  above  solution  can  be  approximated  by 


,♦  =  An  V  -T)  sin  K 


1/3 


T)  sin 


rt 


(-2n)3/2l 


If  K^yg  is  then  replaced  by  the  first  term  of  its  asymptotic  expansion 
(reference  22,  p.  202) 


there  results  finally  for  ^  the  expression 


t  =  B  sin 


fSr  )  X  (-^)"1/4  exp 

4-  (-2 

Vewy 

L  3ew  J 

where  B  is  an  unknown  constant. 


(21) 


The  asymptotic  solution  ( 21 )  makes  it  possible  to  terminate  the 
finite-difference  lattice  at  a  position  on  the  left.  Consider  a  typi¬ 
cal  lattice  point  in  a  column  located  at  rJ=-0  (as,  for  example,  the 
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point  b  in  fig.  l) .  The  neighboring  points 
are  then  as  shown  in  the  adjacent  sketch,  where 
the  point  1+  now  represents  a  ficticious  lattice 
point  located  at  rj"  =  -  ( p+A) .  if  p  is 
taken  sufficiently  large  that  A/p<<l,  then 
it  follows  from  equation  (2l)  that,  to  a  first 
order  of  approximation. 


A _ b\  a 


e„ 


Substitution  of  this  value  of  i|r4  into  the 

previous  equation  (20)  gives  for  the  finite -difference  equation  at 
point  on  the  left-hand  boundary 


1  -  ± 


-  2(l+2p)  *n  =  0 


The  value  of  p  to  be  used  in  any  particular  case  is  determined  on  the 
basiss  of  computational  experience.  In  general,  the  larger  the  value 
of  the  larger  must  be  the  value  of  p  to  assure  that  the  use  of 

the  asymptotic  solution  (2l)  is  justified.  Since  the  over -all  result 
is  insensitive  to  changes  in  the  left-hand  portion  of  the  field,  how¬ 
ever,  the  choice  of  p  is  not  a  critical  matter. 

Points  adjacent  to  the  shock  polar.-  Points  adjacent  to  the  shock 
polar  require  special  treatment  because  of  the  irregularity  of  the 
intervals  encountered  near  the  curved  boundary.  Consider  the  typical 
case  shown  in  the  sketch  (corresponding  to  point  c  in  fig.  l) . 

Here  h  and  k  define  the  length  of  the 
irregular  intervals  relative  to  that  of  the 
regular  interval  A. 


To  obtain  the  desired  accuracy,  it  was 
found  advisable  in  the  present  case  to 
include  three  rather  than  two  neighboring 
points  in  each  of  the  coordinate  directions 
The  value  of  t  at  the  points  2,  4,  and  6 
is  therefore  written 


A _ c  kb./ 

~0/2 
h&  / 


' 3 
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pi  Y~av|  o  . 

~  !0  T)  r[  O  -5*  T)  TJ  TJ  |0 


+  m;  *  |  +  0(a4) 

Q  I  AW  <W  I  '  ' 


to -A  \|r. 


+  A2  * 
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O  *-•  TJ  T] 


A3 

“  “ ™“  Y  Aw  Aw  'V 

O  3  .  T)  T|  -q 


+  0(A4) 


t  n  -  2A  ^ 


♦  (isL*  ♦. 


(2A) 


\ir 

*  «"w  <V  AW 


o  21  q  q  lo  31  qq  q 


+  0(A4) 


These  may  he  looked  upon  as  constituting  three  simultaneous  equations 
for  the  first  three  derivatives  of  t  in  the  horizontal  direction  at 
the  point  0.  Solution  of  these  equations  for  gives 

q  q  ° 


t. 


*np 


l-k 

2+k 


*6 


+ 


2(2-k) 

1+k 


♦4  " 


+  6  + 
k  °  k( 1+k) (2+k)  2J 


+  0(A2) 


is  identical  except  that  k  is 


The  corresponding  expression  for 

yyio 

replaced  by  h  and  +2  ,  ty4,  and  by  ^3,  ti,  and  ^5,  respectively 

Substituting  these  expressions  for  the  two  second  derivatives  into  the 
differential  equation  (10)  and  neglecting  terms  of  0(A2)  then  gives 
for  the  finite -difference  equation  at  0 


k( 1+k) (2+k) 


+  2--2-"—  *4  -  ~  *6  -  2q0  [ 


1+k 


2+k 


[_h(  l+h)(2+h) 


*3  + 


MlrM  t  *  1  _ 

1+h  1  2+h  5J 


3-k 


-  2rf0 


3-h 

h 


to  =  0 


(23) 


This  reduces  to  the  previous  equation  (20)  when  k  =  h  =  1.  (The 
functions  of  h  and  k  which  appear  here  have  been  tabulated  in 
reference  23.  The  intervals  of  tabulation  are  not  always  sufficiently 
small,  however,  to  provide  the  accuracy  needed  in  the  present  work.) 

Points  on  the  shock  polar.-  In  past  applications  of  numerical 
methods  to  problems  involving  curved  boundaries,  it  has  not  ordinarily 
been  the  practice  to  use  a  lattice  with  points  located  on  the  boundary 
itself.  The  prescribed  boundary  conditions  have  then  been  incorporated 
in  the  following  manner  (cf.  references  24  and  25):  First,  the  finite- 
difference  lattice  is  extended,  on  the  basis  of  the  regular  lattice 
spacing,  to  include  ficticious  points  external  to  the  boundary.  This 
makes  the  lattice  geometry  at  internal  points  adjacent  to  the  boundary 
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the  same  as  at  all  regular  internal  points.  Next,  with  the  aid  of  the 
boundary  conditions  and  suitable  interpolation  and  extrapolation  formu¬ 
las,  an  expression  is  obtained  for  the  independent  variable  at  each 
external  point  in  terms  of  the  values  at  neighboring  internal  points. 
Finally,  by  substituting  these  expressions  into  the  finite -difference 
equation  for  a  regular  point,  the  difference  equations  are  written  for 
the  internal  points  adjacent  to  the  boundary.  In  this  way  the  boundary 
conditions  are  incorporated  inplicitly  into  the  difference  equations  at 
internal  points.  The  procedure  is  parallel  in  many  respects  to  that 
used  in  terminating  the  present  lattice  at  the  left-hand  side  of  the 
field. 

Although  a  procedure  of  the  foregoing  type  can  be  devised  to  take 
care  of  the  boundary  conditions  on  the  shock  polar,  a  different  approach 
was  found  advantageous  for  the  present  work.  In  this  approach,  the 
lattice  points  are  placed  directly  on  the  boundary  as  previously 
described,  and  a  difference  equation  is  obtained  at  each  such  point  by 
suitable  finite-difference  approximation  to  the  boundary  condition. 

This  leads  to  a  somewhat  larger  system  of  simultaneous  equations  than 
would  the  more  usual  procedure  but  has  been  found  in  the  long  run  to 
give  more  accurate  results  with  less  total  effort. 

The  boundary  condition  (l4) ,  which  is  thus  the  basis  for  the 
finite-difference  equations  on  the  shock  polar,  can  be  written 


fcw  -  S(tD  tv  =  0 
n  e 


where 


(24) 


The  problem  now  is  to  determine  differ¬ 
ence  expressions  for  the  derivatives 
tv  and  tv  at  points  on  the,  polar. 

11  e 

Consider,  for  example,  the  typical 
situation  shown  in  the  sketch  (corre¬ 
sponding  to  point  d  of  fig.  l) .  To 

determine  the  value  of  \|r  at  each 

0 

of  the  interior  points  1,  2,  and  3  is 
written,  as  before,  in  terms  of  a 
Taylor’s  series  about  the  boundary 
point  0.  Solving  the  resulting  three 
equations  for  the  derivative  1]/^  at  0 
then  gives  0 
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h(2+h)+2(l+h)a  ^  |  ( 1+h) ( 2+h)  t 

h(  1+h)  ( 2+h)  °  2h 


h( 2+h)  ^  +  h(l+h) 
1+h  2  2(  2+h) 


+  0(A3) 


(25) 


This  expression,  which  includes  terms  of  0(a2)>  is  inconsistent  in 
order  of  accuracy  with  the  expressions  previously  employed  in  setting 
up  the  finite-difference  equations  at  internal  points.  Since  end 
differentiation  is,  even  for  a  given  order  of  mathematical  accuracy, 
inherently  less  precise  than  differentiation  at  a  midpoint  (see  error 
terms  in  reference  23) ,  the  retention  of  the  second-order  terms  was  here 
thought  advisable. 


The  determination  of  the  corresponding  expression  for  is  a 

bit  more  involved.  Expanding  at  the  boundary  points  b  and  5  by 
means  of  a  Taylor’s  series  in  two  dimensions  gives 

t4  =  +0-ANU  -  m  *~|  +  ~  +  fAe  ik _ I  +  * _ I  +  0(A3) 

^l0  el°  21  TlT)|o  TJ@  j o  21  eelo 


%  =  %+  Atyvl  +  eA^I  +  %  t — -t  +  eA2* — I  +  +  0(A3) 

llo  0  |o  2.  rjulo  TJ0  | o  2‘  00  lo 

An  expression  for  Is  already  known  from  equation  (25)  in  terms 

0  Jo 

of  ^0,  ty2,  *ond  ^3;  and  an  expression  for  can  similarly  be 

0  0iO 

determined.  The  two  foregoing  expansions  may  thus  be  regarded  as  con¬ 
stituting  two  equations  for  the  three  unknowns  i|r*J  ,  +^^1  ,’  and  ^^1  . 

T)  I  O  T)T)|0  T|  0  j  o 

To  solve  for  ^1  ,  one  more  equation  is  necessary.  This  is  provided 
T)  jO 

by  the  differental  equation  (10),  which  also  applies  on  the  boundary 
and  which  may  be  written  at  the  point  0  as 


i _ 


2~o  =  0 

0  0 


The  solution  for  Itv 

t\  o 


is  then  found  as 
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M  -T  1 

-  -2-  *4 

-5^+0 

+  -L*5 . 

-  —  t^i 

T)|o  A  { 

e+f 

e+f 

e+f 

e+f  0 1, 

1  e-f 


1  /  1  \ 


3+2h 


0(Aa) 


(26) 


where 


e  o 


is  given  by  the  previous  equation  (25).6 


The  required  finite -difference  equation  for  the  point  on  the 
shock  polar  can  now  be  obtained  by  substituting  expressions  (25)  and 
(26)  into  the  boundary  condition  (24)  and  neglecting  the  higher -order 
terms  in  each  case.  The  result  can  finally  be  written 


( 1+h)  ( 2+h)  +2L(  3+2h)  h(  2+h)  +4l(  1+h)  .  h(l+h)+2L(l+2h) 

S  +1 - 1+h - **  + - 2(2+h) - + 


eK  -  fK  xjr, 


h(  2+h)  +2(  1+h) 


6l 


h( 1+h) (2+h)  h( 2+h) 


+  (e-f)K 


^o=0 


(27) 


where  K  =  K  (e,f,tj0)  and  L  =  L  (e,f,rf0)  are  given  by 


(e+f)s(\))+2ef 


1  =  (e-f)  (lo  -  f)  K 

Equation  (27)  is  convenient  for  points  on  the  shock  polar  for 
which  -0.6  ^  T)<0.  For  -1  <*T)0  < -0.6,  the  general  procedure  is  the 


6It  will  be  noted  that  the  coefficients  in  equation  (26)  become 

undesirably  large  as  f  — >  -e  and  are  undefined  when  f  =  -e .  This 
results  from  the  fact  that  the  determinant  of  the  coefficients  in  the 
simultaneous  equations  used  to  obtain  vanishes  when  f  =  -e. 

nl0 

Difficulties  from  this  source  can  be  avoided  by  judicious  choice' of 
the  lattice  points. 
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same  except  that  the  points 
0,  1,  2,  and  3  are  now  more 
conveniently  located  on  a  hori¬ 
zontal  line  and  the  quantities 
e,  f,  and  h  are  redefined 
accordingly  ( see  sketch) .  The 
resulting  finite-difference 
equation  is  identical  with 
equation  (27)  except  that  the 
terms  which  previously  arose 
from  the  expression  for 

6  o 

are  now  multiplied  by  -1. 

Points  on  the  sonic  line.- 


The  difference  equation  for 

points  on  the  sonic  line  is  obtained  by  finite-difference  approximation 
to  the  boundary  condition  (l5)>  wbich  can  be  written,  to  the  accuracy 
required  in  the  numerical  work,  as 


nZ  **  0,0*  ~  ,  , 

iJ0,9)  +  0.342  /  9  -2——-  d0>  =  0  (2 1 

^  J  (e'-e)2/3 

^W 

The  procedure  varies  depending  on  whether  the  difference  equation  is 
being  written  for  the  first  point  below  the  upper  boundary  or  for  one 
of  the  lower  points. 

First  point  below  boundary: 

At  the  first  point  below  the  upper  boundary  the  situation  is  as 
indicated.  To  approximate  the  integral  in  equation  (28),  use  is  made 

of  Guderley's  singular  solution  for  tran¬ 
sonic  flow  over  a  convex  corner  (refer- 

_ B  ence  26) .  On  the  basis  of  this  solution, 

it  can  be  shown  that  the  variation  of  ^ 
along  the  vertical  axis  in  the  immediate 
^  vicinity  of  the  point  B  is  of  the  form 


t  ~  (gL-0)4/3 


If  the  lattice  spacing  is  made  sufficiently  small,  this  asymptotic 
relation  may  be  taken  as  approximately  correct  over  the  entire  interval 
from  B  to  0,  so  that  within  this  interval 
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and 


♦ 


\|r 

0 


4  (M)1/3 

3  *°  A4/3 


For  the  first  point  below  the  boundary,  the  integral  in  equation  (28) 
can  thus  be  written 


+3(0, ^ ~ 

— 2 — - d0’ 

(0'-0)a/3 


^o 

3A4/3 


(Ve')1/3 

[a-(0w-®')]2/3 


^  O  p1  T1/3 

3A£/3  J  (1-t)*/3 


where  T  =  (6w-0' ) /a.  The  integral  on  the  right  can  be  reduced  to 
standard  form  by  means  of  the  substitution  t(1-t)  ==  l/4  z3,  which 
gives  .  * 


tV3 


(1-t)2/Q 


dr 


3  n 1  dz 

^  J0 


This  is  an  elliptic  integral  of  the  first  kind.  Its  value,  as  deter¬ 
mined  from  the  equations  and  tables  of  reference  27,  is 


-1/3 


(1-t) 


2/3 


dT  =  -JL_  x  =  2.650 


,2/3 


1/4 


The  integral  in  equation  (28)  thus  becomes,  in  the  present  case. 


p6  i*(o,0')  „  k  t 

L  ie'  -  3 x  2-650  3* 


0  =  3.533  *° 


2/3 


(30) 


0. 


W 


2h 
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To  approximate  the  derivative  ^  in  equation  (28),  t  is 

n 

expanded  at  points  i  and  ii  by  means  of  a  Taylor's  series  about 

point  Q.  Terms  involving  may  be  omitted  here,  since  the 

T]  Tito 

differential  equation  (lO)  shows  this  derivative  to  be  identically  zero 
at  points  on  the  sonic  line.  The  values  of  t  at  i  and  ii  can 
thus  be  written 


*i  =  *0  -  A 

n 


o  3'-  im 


+  o(a4) 


t±i  =  -  2A 


(2A) 


T1IO  3:  T)T)T)|0 


+  o(a4) 


Solution  of  these  equations  for  ^1  gives,  to  the  second  order  in  A, 

T}  ]  © 


1U 

n|o 


A  U  ° 


(31) 


Substitution  of  expressions  ( 30)  and  ( 31)  into  the  boundary  condi¬ 
tion  (28)  gives  the  following  finite-difference  equation  for  the  first 
point  below  the  upper  boundary: 


*ii  -  /il+  1.208  A1  /3\  io  =  0 

3  6  \6  J 


(32) 


Unlike  the  previous  equations  (20), 
involves  the  value  of  A. 


(22),  (23),  and  (27)  this  equation 
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Lower,  points: 


For  a  general  point  below  the  first  lattice  point,  the  situation 
is  as  represented  in  the  sketch.  The  integral  in  the  boundary  condi¬ 
tion  (28)  is  here  evaluated  in  three  sections.  The 


integral  over  the  lattice  interval  from  B  to  n-1 
is  evaluated  on  the  basis  of  the  asymptotic  relation 
used  before.  The  integral  from  n-1  to  1  is 
evaluated  by  assuming  a  linear  variation  in  \|r  over 
each  of  the  included  intervals  and  then  integrating 
analytically.  The  linear  assumption  is  sufficiently 
accurate  here,  since  this  entire  middle  section  con¬ 
tributes  only  a  relatively  small  portion  of  the 
complete  integral.  The  integral  from  1  to  0, 
which  contains  a  singularity  in  . the  integrand  at 
the  point  0,  is  evaluated  by  expressing  \| f  as  a 
cubic  in  terms  of  its  value  at  the  points  2,  1,  0, 
and  -1  and  then  integrating  analytically  as  before. 
The  added  accuracy  of  the  cubic  is  required  here, 
since  this  last  section  contributes  by  far  the 
majority  of  the  over-all  value.  The  boundary  condi¬ 
tion  (28)  may  thus  be  written 

(0,3)  +  0.342  ( Ji+J2+J3)  =  0  (33) 

where  the  J*  s  represent  the  three  component 
integrals  just  described. 


Proceeding  to  the  details  of  the  above  procedure,  the  integral 
from  B  to  n-1  is  first  written 


,  .  f* >  *1<°.3'>  A.  Wn-1  r1  Tl/3  __ 

1  4  InA-tV?')!*'3  "3A l  (n-T)*/3 

yw 

where,  as  before,  t  =  (0w-0*)/ A.  This  integral  can  be  expressed,  if 
desired,  as  the  difference  of  two  elliptic  integrals  of  the  first  kind. 
For  present  purposes,  however,  it  is  more  convenient  (and  sufficiently 
accurate)  to  expand  (n-T)~2/3  according  to  the  binomial  theorem  and 
integrate  termwise.  This  gives  finally 


Ji  = 


(nA)  2/3 


1— 

+  —  +  —  + 

21n  9ne 


160  55 

-  +  -  + 

1053n3  W36n4 


This  expression  is  used,  of  course,  only  for  n^2. 
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On  the  basis  of  the  assumption  of  a  linear  variation  of  ♦  between 
adjacent  lattice  points,  the  integral  from  n-1  to  1  becomes 


t  =  rew-(n-l)A  tfr(0>) 

2  [nA-(ew-0’)]£/3 


'V 

d0'  = 


eV-A 


W+m  P9' w-(n-m)A 


A  J*'  ,  ,  [nA-(0w-0')]2/a 

m=l  0w-(n-m-l)A 


d0 


Carrying  out  the  integration  gives 

n-2 

J2=^  £  [(m+D1/3-ml/3  ]  (v+m+i)  (35a) 

m=i 

This  expression  is  valid  for  n  (For  n  =  2,  J2  obviously  does 

not  exist.)  For  n  >3  it  is  convenient  to  rewrite  the  summation  so 
that  the  value  of  ♦  at  a  given  point  is  not  repeated  in  successive 
terms  of  the  series.  This  is  done  by  separating  expression  (35a)  into 
two  series  (one  with  and  one  with  +in+i)  >  expanding  these  series, 

and  then  regrouping  terms.  The  result  is  finally 


21/3-! 


h 


+ 


n-2 


l 


(m+l)  1/3-2m1/3  +  (m-l) 1/3 


+m' 


(n-l)1/a-(n-2)1/3  tn.xJ 


(35b) 


This  expression  is  valid  for  n  ^ 

To  evaluate  the  integral  from  1  to  0,  ♦  is  represented  within 
this  interval  by  a  cubic  of  the  form 
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*s 


1 


* 


where  a,  b,  c,  and  d  are  determined  such  that  i|r  has  the  proper 
values  at  the  points  2,  1,  0,  and  -1.  This  expression  is  to  be  substi¬ 
tuted  into  the  integral 


w 


•nA  ^(0,0') 

[nA-iev-e'W7* 


ae' 


The  result  is  finally,  after  evaluation  of  the  coefficients  a,  b,  c, 
and  d, 


j3 


-1 

A2/3 


*  +  13 
28  1  lA 


*o  + 


(36) 


The  finite-difference  equation  for  a  general  point  on  the  sonic 
line  can  no w  be  obtained  by  replacing  the  JTs  in  equation  (33)  hy  the 
expressions  ( 3*0  9  (35)*  and  (36)  and  using  the  previous  expression  (31) 
for  the  derivative  The  result  is  a  lengthy  linear  equation 

T)  I  o 

involving  the  value  of  \| f  at  the  points  i,  ii,  -1,  0,  1,  2,  . .., 
n-1,  n.  Fortunately  for  the  later  relaxation  work,  the  coefficients 
of  the  terms  turn  out  to  be  relatively  small  for  all  points  above  the 
point  2. 


Distribution  of  mesh  points.-  When  an  attempt  is  made  to  solve  the 
present  problem  with  a  coarse  m'esh,  it  is  soon  found  that  most  of  the 
variation  in  ^  takes  place  in  a  relatively  small  region  near  the 
intersection  of  the  shock  polar  and  sonic  line.  To  obtain  a  suffi¬ 
ciently  accurate  solution  in  a  practicable  length  of  time,  it  is 
therefore  necessary  to  employ  a  graded  lattice,  that  is,  a  lattice 
which  has  different  spacing  in  different  parts  of  the  field.  Figure  2 
shows  the  distribution  of  lattice  spacing  found  satisfactory  in  a 
typical  case  ( =  1.6).  The  particular  arrangement  shown  here 
involves  a  total  of  228  lattice  points.  For  other  values  of  *0W,  the 
grading  of  the  lattice  follows  the  same  general  scheme.  Obviously, 
however,  the  total  number  of  lattice  points  must  be  increased  as  the 
upper  boundary  is  moved  farther  from  the  shock  polar. T 

Formulas  (3*0,  ( 35)  >  and  (36),  which  are  used  to  approximate  the 
integral  along  the  sonic  line,  presume  the  existence  of  lattice  points 
at  a  uniform  interval  over  the  full  distance  from  the  upper  boundary 
to  the  point  in  question.  This  condition  is  not  fulfilled  in  a  graded 


^Occasionally,  when  two  points  on  or  adjacent  to  the  shock  polar  fall 
very  close  together,  one  of  the  points  is  arbitrarily  omitted.  An 
expression  for  the  omitted  value  of  which  is  then  necessary  to 
complete  the  difference  equation  at  neighboring  points,  is  found  by 
parabolic  interpolation  between  the  values  at  the  available  locations • 


\ 
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lattice  such  as  that  indicated  in  figure  2,  Some  modification  of  the 
method  must  therefore  be  made  to  obtain  the  finite-difference  equation 
for  a  point  on  the  sonic  line  in  one  of  the  regions  of  finer  mesh. 

This  requirement  was  satisfied  by  means  of  a  simple  averaging  process 
in  which  the  contribution  of  nonexistent  fine-mesh  points  is  replaced 
by  an  average  contribution  expressed  in  terms  of  t  at  bracketing 
points  on  the  available  coarser  net.  Since  the  contribution  of  indi¬ 
vidual  points  is  small  even  for  points  only  moderately  removed  from 
that  at  which  the  equation  applies,  a  rather  crude  averaging  process 
is  sufficient  in  most  cases.  (The  details  need  not  be  given  here  as 
they  would  soon  become  apparent  to  anyone  working  with  the  method.) 

When  the  averaging  procedure  would  not  be  sufficiently  accurate  (as 
when  the  point  at  which  the  equation  applies  is  near  the  line  of  demar¬ 
cation  between  two  different  sized  meshes),  fictitious  intermediate 
points  are  introduced  into  the  coarser  net  and  the  value  of  ^  at 
these  points  is  obtained  from  plots  of  the  distribution  of  f  along  the 
sonic  line. 


Solution  of  Finite-Difference  Equations 


By  the  methods  of  the  foregoing  section,  a  finite-difference 
equation  can  be  obtained  for  each  lattice  point  in  the  hodograph  plane. 
The  result  is  a  large  number  of  simultaneous  algebraic  equations 
involving  an  equal  number  of  unknown  values  of  i|r.  Since  the  number 
of  unknowns  in  each  equation  is  small,  the  equations  lend  themselves 
well  to  solution  by  relaxation  techniques.8 

The  mechanics  of  the  relaxation  process  have  been  well  described 
by  various  authors  (references  15,  16,  17,  and  18)  and  need  not  be  gone 
into  here.  For  present  purposes  it  was  found  satisfactory  to  take 
in  the  boundary  condition  (l6)  equal  to  10,000  and  work  with  integer 
values  of  if  throughout  most  of  the  field.  The  residuals  in  the  relax¬ 
ation  process  were  eliminated  to  within  limits  of  ±2  (with  due  care, 
of  course,  that  all  residuals  in  any  given  area  were  not  predominately 
of  the  same  sign) .  To  obtain  satisfactory  smoothness  of  the  solution 
near  the  left-hand  boundary  in  some  examples,  it  was  necessary  in  this 
region  to  work  with  values  of  \|r  to  0.1  and  eliminate  residuals  to 
within  ±0.5.  Whenever  the  coefficients  in  the  finite-difference 


8 It  is  interesting  to  note  that,  of  the  complete  set  of  simultaneous 
equations,  only  two  —  those  for  the  points  on  the  shock  polar  and 
sonic  line  immediately  adjacent  to  the  point  E  —  are  not  homo¬ 
geneous.  Only  this  fact  prevents  the  solution  of  the  complete  set 
from  being  identically  zero. 
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equations  were  relatively  small,  the  corresponding  terms  were  neglected 
in  the  point -by -point  adjustment  of  i|f .  The  error  so  introduced  was 
eliminated  periodically  by  recomputing  the  residuals  using  all  terms  in 
the  finite-difference  equations.  This  procedure  was  particularly  help¬ 
ful  in  the  case  of  the  lengthy  equations  which  apply  at  points  on  the 
sonic  line.  The  transition  between  the  various  regions  of  the  graded 
lattice,  which  is  not  often  discussed  in  the  literature,  was  accom¬ 
plished  by  the  use  of  overlapping  fields  in  essentially  the  manner 
described  in  reference  28. 

By  means  of  the  foregoing  procedures,  the  boundary- value  problem 
in  the  holograph  plane  has  been  solved  for  values  of  $w  of  1.3,  1.6, 

2. It-,  and  4.2.  These  values  correspond,  respectively  (see  equation  (l8)), 
to  values  of  the  transonic  similarity  parameter  |0  of  I.O58,  0.921, 

0* 703*  and  0.484  as  given  previously  in  reference  1.  As  shown  in  the 
earlier  paper,  calculations  for  these  values  were  sufficient  to  bridge 
the  gap  between  the  previous  results  of  Guderley  and  Yoshihara 
at  £0  =  0  (reference  29)  and  the  analytic  results  which  can  be  obtained 
when  the  bow  wave  is  attached  and  the  flow  is  completely  supersonic 
(reference  30).  As  an  example  of  the  solution  in  the  hodograph  plane, 
the  variation  of  for  Qv  =  1.6  is  shown  as  a  function  of  r\  and 
0  in  figure  3.  (These  results  correspond  to  the  results  shown  in  the 
physical  plane  in  fig.  2  of  reference  1.)  This  figure  shows  clearly 
the  rapid  variation  of  near  the  intersection  of  the  shock  polar 
and  sonic  line.  The  calculated  values  of  corresponding  to  figure  3 
are  listed  in  table  I  at  the  end  of  the  report. 


TRANSFORMATION  TO  PHYSICAL  PLANE 
Flow  Field 


The  transformation  from  the  hodograph  plane  to  the  physical  plane 
is  governed,  in  the  small-disturbance  theory,  by  the  following  equa¬ 
tions  (cf.  reference  5): 


dx  = 


-~+1^  (n  ie  dt)  +  de) 

pa  1 


dy 


l 

P*a* 


(^  d-n  +  i|r0  d0) 


-=—  dt 
p*a* 


The  second  of  these  equations  implies  that,  in  a  flow  field  determined 
according  to  the  transonic  small-disturbance  theory,  all  streamlines 
appear  as  straight  lines  parallel  to  the  horizontal  axis.  When 
expressed  in  terms  of  t)  and  9  the  foregoing  equations  become 
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dx  =  (2Tf  U  dr,  +  U  d Q) 

P*a*  v  2  e  n 

dy  =  — drf  +  ^  d?)  =  — —  d  \|r 
P*a*  T)  0  P*a* 

The  length  Z  of  the  wedge,  which  is  equal  to  one-half  the  chord 
of  the  double-wedge  profile,  can  be  found  by  integrating  the  first  of 
these  equations  over  the  upper  boundary  OB  of  the  hodograph.  (See 
sketch.)  This  gives 


0 


With  this  relation  the  previous  transformation  equations  can  he  put  in 
the  dimensionless  form 

d  (-)=  (29f  ^  dSf  +  ^  d0)  (37) 

\cy  «v  9  n 

( 7+1) 1/3  d(k)  =  —  di  (38) 

V  2  \c  /  4lw 
where  Iw  represents  the  integral 

Iw  =  /  T  (ti^w)  dT) 

J  e 


— co 


(39) 
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To  obtain  the  flow  field  in  the  form  given  in  reference  1,  equation  ( 38) 
must  be  rewritten  in  terms  of  the  ordinate  function 


1/3  /y 


Y  =  [(r+l)(t/c) ] 


The  result,  derived  with  the  aid  of  the  relation 


n02/3 


d?  a  [(r+i)(t/c)J1/3  iQf)  - 


(29J1/3 


( 40a) 


Integrating  this  relation,  subject  to  the  condition  that  Y  =  0  when 
t  =  0,  gives 


Y  =  [(7+l)(t/c)]1/3 


(29v)1/3 


(1+Ob) 


To  utilize  the  foregoing  equations  for  actual  computations,  it  is 
first  necessary  to  evaluate  Iw.  Since  numerical  values  of  \|r  are 
available  in  the  holograph  only  for  -p=T]^0,  the  evaluation  must  be 
carried  out  in  two  parts  as  follows: 


n  W,K)  a?  +  T~* 


0 


The  first  integral  is  evaluated  from  the  results  of  the  numerical 
solution  by  mechanical  integration  of  a  curve  of  r\  ^(t),0w)  versus  *rj 

6 

The  values  of  the  derivative  used  for  this  purpose  are  obtained  from 
the  equation 


'M  ti  +4^2 
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where  the  notation  is  as  shown  in  the  sketch. 

This  equation  is  derived  in  the  same  way  as  equa¬ 
tion  (3l)>  except  that  +0  is  here  taken  equal 
to  zero  in  accord  with  the  boundary  condition.9 
It  can  be  shown  from  Guderley's  singular  solution 
for  corner  flow  that  for  small  negative  values 

of  ?J  the  curve  of  ?f  i|r#>*(q,0w)  must  behave  essen- 

Q 

tially  as  |rf|3/2.  This  result  is  useful  in 
fairing  the  numerical  results  near  q  =  0.  The 
first  integral  in  equation  (4l)  contributes  by 
far  the  majority  (about  99  percent)  of  the  total 
value  of  4, 


To  evaluate  the  second  integral  in  equation  (4l)  use  in  again  made 
of  the  asymptotic  solution  (21).  For  this  purpose,  the  constant  B  is 
determined  such  that  the  value  of  'L,  given  by  the  asymptotic  solution 

e 

matches  the  numerically  determined  value  at  the  point  (-|3,@w). 
Substitution  of  equation  (21)  into  the  second  integral  of  equation  (4l) 
then  gives 


r  "|3~  p 1/4  *&;-(>  A) 

/  1  WllA)  dll  = - - - 

J  __  6  r 


-  -Sr  (2P)3/2] 
L  3©w  J 


exp 

.-P 


(_rf)3/4  exp 


-  00 


-  <-2tD3/2l  Hi 

L  3^w  J 


(**■3) 


where  1hv(-(3,0w)  determined  from  equation  (42)  applied  at  q  =  -(3. 
0 

The  integral  on  the  right  is  transformed  through  the  substitution 


(-21I)3/2  =  » 

3f>w 


^he  fact  that  the  second  derivative  may  be  taken  as  zero  in 

the  present  derivation  follows  from  the  boundary  condition  and  the 
differential  equation  (l0)« _ _ 
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which  gives 
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f*  < 
J  -  00 

_~)3/4  eXp 

JL  (-2n)3/2!  <fli  = 

.  3@w  J 

1 

/30w  y/e 

* 

•(Jji/6e"wd m 

-4-  (2s)3'2 

3©w 

3x23/4 

\TJ  o 

The  integral  he're  is  then  found  with  sufficient  accuracy  by  means  of 
the  following  asymptotic  formula,  valid  for  large  values  of  the  lower 
limit  (see  pp.  95-96  of  reference  31): 

CD 


V-l  e“u  aw  =  z^1  e~z 


3 


Equation  (43)  thus  reduces  finally  to 


/ K  4?-  -  Q  J/a  *~(-<s,3»)  (H) 

With  the  value  of  Iw  known,  equations  (37)  and  (ko)  can  be  used 
to  obtain  the  coordinates  x/c  ^and  Y  corresponding  to  any  point  in 
the  T), 9  plane.  The  value  of  Y  is  obtained  by  direct  substitution 
of  the  appropriate  value  of  +  into  equation  (40b).  .  The  value  of  x/c 
must  be  found  by  suitable  integration  of  equation  (37).  The  location 
of  the  vertex  A  of  the  shock  wave  is  found,  for  example,  by  integrat¬ 
ing  equation  (37)  along  the  line  OA  in  the  holograph.  (See  sketch 
on  p.  30.)  If  the  leading  edge  0  is  taken  as  the  origin  in  the 
physical  plane,  this  gives* 

(t  \  m  f  ^  V^0’ m 

-  00 

The  integral  here  is  evaluated  in  two  parts  following  the  procedure 
previously  used  in  determining  Iw. 

For  the  abscissa  of  a  point  F  on  the  shock  wave,  equation  (37) 

gives 

=  (l)  +  A-  f  (2?i  ^  dn  +  ae) 

\7A  «v  J.  e  ti 
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where  the  integration  is  now  taken  from  A  to  F  along  the  shock 
polar.  For  purposes  of  numerical  evaluation,  the  integrand  here  can  be 
simplified  by  writing 


2t)  i|r^  dr)  +  ^ 


9  il 


d@ 


dq  +  K  d0 

11  e 


d\|r 


=  21+  (trfA^)(de/d^)  ^ 

(t^/iL)  +  (d?/aii) 

T]  9 

If  and  (dQ/dq)  are  replaced  by  the  appropriate  functions  of  rj 

tj  e 

from  equations  (l4),  there  results  finally 


The  integral  in  this  equation  is  evaluated  by  plotting  a  curve  of  \|r 
versus  from  the  numerical  results  along  the  shock  polar  and 

carrying  out  the  necessary  integration  by  mechanical  means. 


The  abscissa  of  a  point  on  the  sonic  line  is  found  by  integrating 
equation  (37)  along  the  6  axis  from  B  to  C.  Since  point  B  is 
situated  in  the  physical  plane  at  x/c  =  l/2,  this  gives 


1  _1_ 

2  +  4lw 


1M0/8)  de 
il 


(^7) 


The  integral  here^ is  evaluated  by  mechanical  integration  of  a  curve  of 
versus  9,  where  t^(O,0)  is  found  from  equation  (31).  As  can 

T|  T)  ^ 

be  seen  from  equation  (28)  and  relation  (29),  ^(0,0)  in  the  vicinity 

of  point  B  varies  essentially  as  (£w^)2^3-  This  fact  is  of  use  in 
drawing  the  curve  of  ^(0,0)  near  0=9W.  It  can  further  be  seen  with 

the  aid  of  equation  (40  b)  that  near  the  shoulder  of  the  wedge  the  trans¬ 
formed  sonic  line  has  the  form 


This  relation  is  useful  in  establishing  the  detailed  shape  of  the  sonic 
line  in  the  physical  plane.  It  shows,  in  particular,  that  the  sonic 
line  will  have  a  vertical  tangent  and  an  infinite  curvature  at  the 
shoulder  of  the  wedge. 
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Pressure  Distribution  and  Drag 


To  complete  the  analysis  of  the  front  wedge,  it  is  left  to  deter¬ 
mine  the  pressure  distribution  and  drag.  Integration  of  equation  (37) 
gives  for  the  chordwise  location  on  the  wedge  of  a  given  value  of  rf 


or 


x  _  1 

c  2Iw 


i  a? 

e 


m 


The  speed  parameter  |  =  (M2-l) /[(y+l) ( t/c) ]2/3  which  was  used  to 
present  the  results  in  reference  1,  is  related  to  r)  by  the  following 
equation,  derived  from  equations  (20)  and  (21)  of  reference  1  and 
equations  (l),  (9),  and  (l8)  of  the  present  paper: 


6  =  ^o  *  =  ~ 


2i/s 


q  2/3 

yw 


(^9) 


With  these  equations,  the  distribution  of  i  as  a  function  of  x/c 
is  readily  determined.  The  integration  of  equation  (48)  is  carried 
out  by  mechanical  means  using  the  same  curve  previously  employed  to 
determine  Iw.  To  fair  the  resulting  |  curve  in  the  vicinity  of  the 
shoulder,  use  is  again  made  of  Guderley* s  analytical  findings,  which 
show  that  in  this  vicinity 


The  chordwise  distribution  of  the  generalized  pressure  coefficient 


(7+l)1/3/(t/c)2/3 


is  related  to  the  distribution  of  £  by  the  following  equation  taken 
from  pages  l6  and  23  of  reference  1: 

c p  =  -  2  (i-i0) 


(50) 


36 


HACA  TIT  2588 


To  find  the  portion  of  the  generalized  drag  coefficient 
~d  =  [(7+l)/(t/c)5/3]  cd 

contributed  by  the  front  half  of  the  double-wedge  profile,  equation  (13a) 
of  reference  1  gives 


(51) 


> 


To  allow  for  the  singularity  in  Cp  at  the  leading  edge,  the  integra¬ 
tion  is  best  carried  out  directly  in  terms  of  the  holograph  variables. 
To  this  end,  substitution  from  equations  (37),  (^9),  and  (50)  gives 


~2 

T1 


f  'V*  .  r\t 

t^(T)>6w)  dTJ 

e 


1 

j 


(52) 


The  integration  here  is  performed  in  two  parts  following  procedures 
parallel  to  those  used  in  evaluating  Iw 


CHARACTERISTICS  CONSTRUCTION  OVER 
REAR  OF  AIRFOIL 


The  characteristics  in  the  plane  (t]>0)  are  given  by  the 

following  relation  obtained  from  equations  (3)  and  (9): 

Q  =  const.  ±  - -  rf  3//2  (53) 

3 

The  corresponding  directions  of  the  Mach  lines  in  the  generalized 
physical  plane,  as  determined  from  this  relation  and  the  transformation 
equations  (37)  and  (40a),are 


dY 


3f  ^/3 


=  + 


(29w) 


d(x/c)  ( 2rj) 1//2 


(5M 
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To  the  present  order  of  approximation,  therefore,  the  slope  of  the  Mach 
lines  is  independent  of  the  local  inclination  0.  '  This  is  consistent 
with  the  previous  result  that  the  streamlines  must  appear  in  the  physi¬ 
cal  plane  as  horizontal  straight  lines.  As  a  consequence,  the  construc¬ 
tion  of  the  Mach  net  over  the  rear  of  the  airfoil  is  particularly 
simple  in  the  small-disturbance  theory. 

To^aid  in  the  construction,  the  equations  for  the  characteristics 
in  the-  T),0  plane  can  be  conveniently  written  in  the  form 

0=6^  +  2^3/2 

3 

0  =  gd  -  9J3/2 

3 


The  symbols  ©u  and  ©d  denote,  respectively,  the  ordinates  at  which 
the  upgoing  and  downgoing  characteristics  through  a  given  point  (r),0) 
intersect  the  vertical  axis.  Elimination  of  0  between  these  equa¬ 
tions  gives 


n  = 


32/3 

- -  (©  d 

25/3  V  d 


eu) 


2/3 


which  can  be  substituted  into  equation  (54)  to  obtain 


dY  =  +  /4gw/3\  i/3 

d(x/c)  ~  \  ®dr®u/ 


(55) 


(56) 


This  is  the  basic  relation  for  the  characteristics  construction  in  the 
physical  plane. 

The  construction  of  the  Mach  net  itself  follows  a  simple  lattice- 
point  procedure  (cf.  reference  32)^  By  identifying  each  Mach  line  with 
its  appropriate  value  of  0U  or  ©d,  the  value  of  dY/d(x/c)  at  the 
intersection  of  any  two  Mach  lines  can  easily  be  determined  from 
equation  (56)  (or  its  graphical  equivalent).  The  basic  construction 
necessary  to  locate  an  unknown  point  c  from  the  location  of  two  known 


O^CD 
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points  a  and  b  is  then  as  indicated  in  the  sketch.  The  construction 
proceeds  rapidly  since,  as  pointed  out,  variations  in  the  inclination  of 


flow  need  not  be  considered  in  establishing  the  direction  of  the  Mach 
lines .  Where  desired,  the  value  of  6  can  be  found  from  the  relation 


+  0 


u 


(57) 


The  corresponding  value  of  %  is  given  by  equation  (55). 10 

Figure  4  shows  a  typical  Mach  net  constructed  by  stepwise  appli- 
ation  of  the  foregoing  procedure.  This  net  is  for  the  case  of 
w  =  1*6  (lo  =  0.921)  and  corresponds  to  the  flow  field  shown  for  the 
front  of  the  airfoil  in  figure  2  of  reference  1.  The  construction  is 


10In  practice,  the  ^construction  is  actually  carried  out  most  easily  in 
a  plane  of  Y/(20w)1/3  versus  x/c  with  the  slope  of  the  Mach 
lines  given  by 

d[Y/(2?w)1/a]  _  +  /  2/3  \1/3 
a(x/c)  \ed-©u/ 

This  allows  a  single  graph  of  slope  versus  (Sd-Su)  to  suffice  for 
all  values  of  It  also  provides  somewhat  more  convenient  pro¬ 

portions  for  the  construction  of  the  Mach  net. 
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begun  at  the  shoulder  of  the  airfoil  (x/c  =  1/2,  Y  =  0)  with  the  values 
of  0 u  selected  to  provide  approximately  equal  spacing  between  the  Mach 
lines  of  the  expansion  fan.  From  the  shoulder,  the  construction  is 
carried  outward  to  the  sonic  line  and  then  inward  to  the  rear  surface 
of  the  airfoil.  The  drawing  of  the  Mach  line  segments  adjoining  the 
sonic  line  might  appear  at  first  to  offer  some  difficulty,  since  a 
linear  average  is  obviously  unreliable  to  determine  a  mean  inclination 
in  this  vicinity.  Actually,  no  trouble  is  encountered  from  this  source, 
since  the  point  at  which  each  Mach  line  meets  the  sonic  line  is  already 
known  from  the  hodograph  solution  for  the  subsonic  field.  The  construc¬ 
tion  of  the  last  segment  approaching  the  sonic  line  thus  reduces  to  a 
matter  of  simply  connecting  two  known  points.  The  slope  of  the  first 
segment  leaving  the  sonic  line  is  found  by  either  (a)  multiplying  the 
slope  of  the  approaching  segment  by  -1,  or  (b)  determining  a  mean 
inclination  based  on  the  easily  demonstrated  fact  that  a  Mach  line  in 
the  vicinity  of  the  sonic  line  behaves  essentially  as  a  semicubical 
parabola.11  It  is  immaterial  to  the  final  result  which  procedure  is 
used.  The  identity  of  the  Mach. lines  reflected  from  the  rear  surface 
of  the  airfoil  is  determined  from  equation  (57)  plus  the  boundary  condi¬ 
tion  that  at  this  surface  £  =  As  can  be  seen  by  comparing 

figure  .4  of  the  present  report  with  figure  2  of  reference  1,  only  a 
relatively  small  portion  of  the  sonic  line  need  be  known  to  determine 
conditions  on  the  rear  of  the  airfoil. 


REMARKS  ON  ACCURACY  OF  SOLUTION  * 


Quantitative  statements  with  regard  to  the  accuracy  of  the  present 
results  are  difficult  to  make.  Fortunately,  however,  a  check  on  the 
accuracy  of  the  solution  is  available  in  the  work  itself.  This  check 
derives  from  the  fact  that,  in  the  subsonic  portion  of  the  field,  the 
calculated  location  of  a  given  velocity  in  the  physical  plane  should, 
theoretically,  be  independent  of  the  path  of  integration  which  is 
followed  in  the  hodograph.  Thus,  for  example,  the  position  of  the 
velocity  *rj  =  0,  *8  =  1,  which  defines  the  point  of  intersection  E  of 
the  shock  wave  and  sonic  line,  should  be  the  same  irrespective  of 
whether  it  is  found  from  equations  (45)  and  (46) 

(58) 


11 The  latter  possibility  was  pointed  out  to  the  authors  by 
Gottfried  Guderley. 


-1 


E 


=  —  rf  f  (rf,0)  dr(  + 

^1w  .  /  n 


t 


E 
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1+T)  d* 


w 


or  from  equation  (47) 
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=  i + 77-  r  &  (59) 

2  4IW  J  T) 

ew 

Actually,  as  was  observed  in  connection  with  figure  1  of  reference  1, 
the  results  of  the  two  determinations  show  a  small  discrepancy.  Such  a 
result  would  be  expected  in  any  finite-difference  solution. 

Discrepancies  of  the  type  noted  can  arise  from  two  sources:  (l) 
numerical  inaccuracies  in  the  relaxation  solution  of  the  finite- 
difference  equations  or  in  the  transformation  to  the  physical  plane;  (2) 
inaccuracies  caused  by  the  fact  that  the  finite-difference  equations 
themselves  are  not  an  exact  representation  of  the  boundary- value  problem 
for  the  original  partial  differential  equation.  Experience  with  various 
refinements  in  the  calculations  indicates  that  the  discrepancies  here 
are  primarily  of  the  latter  origin.  Early  computations  with  a  coarse 
lattice  and  relatively  crude  finite-difference  equations  showed  a  con¬ 
siderable  gap  between  the  shock  wave  and  the  end  of  the  sonic  line. 
Increasing  refinements  in  the  grading  of  the  lattice  and  in  the  deriva¬ 
tion  of  certain  of  the  finite-difference  equations  gave  progressive 
improvement  in  reducing  this  gap.  This  improvement  came  about  primarily 
as  a  result  of  progressive  reduction  in  the  value  of  the  integral  Iw, 
the  other  integrals  in  equations  (58)  and  (59)  being  relatively  unaf¬ 
fected  by  the  refinements  in  the  calculations.  Indications  are  that, 
in  the  results  which  were  taken  as  final,  the  values  of  I^A Tl,^w)  I 

Q 

and  hence  of  Iw  are  still  somewhat  too  large.  This  means  (see 
equation  (40b))  that  the  ordinates  of  the  shock  wave  and  sonic  line  are 
probably  somewhat  smaller  than  they  should  be.  The  same  is  probably 
true,  in  general,  of  the  corresponding  values  of  | x/c I .  Calculations 
of  the  chordwise  distribution  of  tj  on  the  surface  of  the  airfoil  are, 
however,  considerably  more  precise,  since  the  errors  in  the  two  inte¬ 
grals  in  equation  (48)  tend  to  compensate.  The  refinements  in  the 
computations  were,  in  fact,  carried  to  the  point  where  further  better¬ 
ment  caused  only  negligible  change  in  the  pressure  distribution  and 
over-all  drag.  Further  evidence  of  the  accuracy  of  the  results  in  this 
regard  is  provided  by  the  ease  with  which  the  computed  values  fair  into 
the  results  of  Guderley  and  Yoshihara  at  =  0  and  into  the  analyti¬ 
cal  curves  which  are  available  when  the  bow  wave  is  attached  and  the 
flow  is  completely  supersonic.  (See  figs.  4  and  6  of  reference  1.) 
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TABLE  I.-  CALCULATED  VALUES  OF 


=  1.6  (|0  =  0.921) 
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Figure  l-  Illustrative  finite -difference  lattice  in  the  plane. 


Figure  2.  -  Distribution  of  lattice  spacing  for 
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